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Abstract. Classical ellipsoidal and sphero-conal harmonics are polynomial solutions of 
the Laplace equation that can be expressed in terms of Lame polynomials. Generalized 
ellipsoidal and sphero-conal harmonics are polynomial solutions of the more general Dunkl 
equation that can be expressed in terms of Stieltjes polynomials. Nivcn's formula connec- 
ting ellipsoidal and sphero-conal harmonics is generalized. Moreover, generalized ellipsoidal 
harmonics are applied to solve the Dirichlet problem for Dunkl's equation on ellipsoids. 
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1 Introduction 

The theory of ellipsoidal and sphero-conal harmonics is a beautiful achievement of classical 
mathematics. It is not by accident that the well-known treatise "A Course in Modern Analysis" 
by Whittaker and Watson (20] culminates in the final chapter "Ellipsoidal Harmonics and Lame's 
Equation" . 

An ellipsoidal harmonic is a polynomial u(xq, xi, . . . , x^) in k+l variables xq, x\, . . . , x k which 
satisfies the Laplace equation 

A-=Eg = (1.1) 

and assumes the product form 

u(x , Xl , ...,x k ) = E(t Q )E{t{) ■ ■ ■ E(t k ) (1.2) 

in ellipsoidal coordinates (in,ii, . . . ,tk) with E denoting a Lame quasi-polynomial. 

A sphero-conal harmonic is a polynomial u(xq,x\, . . . ,Xk) which satisfies the Laplace equa- 
tion (jLlj) and assumes the product form 

u(x , x u ...,x k ) = r m E{ Sl )E{s 2 ) ■ ■ ■ E{s k ) (1.3) 

in sphero-conal coordinates (r, s±, S2, ■ ■ ■ , Sfe)- Again, E is a Lame quasi-polynomial. 

In most of the literature, for example, in the books by Hobson |S] and Whittaker and Wat- 
son |20j . ellipsoidal and sphero-conal harmonics are treated as polynomials of only three vari- 
ables. The generalization to any number of variables is straight-forward. Since we plan to work 

*This paper is a contribution to the Vadim Kuznetsov Memorial Issue "Integrable Systems and Related Topics". 
The full collection is available at http://www.emis.de/journals/SIGMA/kuznetsov.html 
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in arbitrary dimension we will employ ellipsoidal and sphero-conal coordinates in algebraic form; 
see Sections 3 and 4. In IR 3 we may uniformize these coordinates by using Jacobian elliptic func- 
tions. Sphero-conal harmonics are special instances of spherical harmonics. Spherical harmonics 
in any dimension can be found in books by Hochstadt [7j and Miiller |15j . 

We will extend the theory of ellipsoidal and sphero-conal harmonics by replacing the Laplace 
equation by the equation 

k 

&gU ■= jV'jU = 0, a = (a ,ati,...,a k ) (1.4) 

3=0 

introduced by Dunkl In (|1.4|) we use the generalized partial derivatives 

d u(x) — u(ajx) 
VjU(x) := ^— u{x) + ctj — , x = (x ,xi,...,x k ), (1.5) 

OXj Xj 

where aj is the reflection at the jth coordinate plane: 

aj(x ,xi, ... ,x k ) = (x ,xi, . . . ,Xj-i, -xj,x j+ i, . . . ,x k ). 

Equation (|1.4[) contains real parameters oto,ai, . . . ,a k . If aj = for all j then the equation 
reduces to the Laplace equation. 

A generalized ellipsoidal harmonic is a polynomial u(xq,xi, . . . ,x k ) which satisfies Dunkl's 
equation (|1.4j) and assumes the product form (|1.2|) in ellipsoidal coordinates but with E now 
denoting a Stieltjes quasi-polynomial. Generalized ellipsoidal harmonics will be treated in Sec- 
tion 3 while Stieltjes quasi-polynomials are introduced in Section 2. 

A generalized sphero-conal harmonic is a polynomial u(xq, x%, . . . , x k ) which satisfies Dunkl's 
equation (|1.4|) and assumes the product form (|1.3|) in sphero-conal coordinates. Again, E is 
a Stieltjes quasi-polynomial. Generalized sphero-conal harmonics will be considered in Sec- 
tion 4. 

It is very pleasing to see Stieltjes quasi-polynomials taking over the role of Lame quasi- 
polynomials. Stieltjes polynomials have been considered for a long time but they did not appear 
in the context of separated solutions of the Laplace equation. Therefore, our paper shows how 
Stieltjes polynomials become part of the theory of "Special Functions". 

It is quite remarkable that all the known results for classical ellipsoidal and sphero-conal 
harmonics carry over to their generalizations. In Section 5 we generalize formulas due to Hobson 
[13 Chapter 4]. In Section 6, as a consequence, we prove a generalization of Niven's formula 
|2(Jl Chapter 23] connecting ellipsoidal and sphero-conal harmonics. In Section 7 we apply 
generalized ellipsoidal harmonics in order to solve a Dirichlet problem for (jl.4[) on ellipsoids. 
This generalizes the classical result that ellipsoidal harmonics may be used to find the harmonic 
function which has prescribed values on the boundary of an ellipsoid. Finally, we give some 
examples in Section 8. 

We point out that parts of Section 4 overlap with the author's paper ^j]- The contents of the 
present paper are also related to the book by Dunkl and Xu and the papers by Liamba |14j 
and Xu |21j . however, these works do not involve Stieltjes polynomials. We also refer to papers 
by Kalnins and Miller [HJ 03 US]. The paper ^0] addresses Niven's formula from a different 
perspective. Kutznetsov ^H] an d Kutznetsov and Komarov ^2] have also worked in related 
areas. Kutznetsov jointly with Sleeman wrote the chapter on Heun functions for the Digital 
Library of Mathematical Functions. Stieltjes polynomials appear in this chapter. 
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2 Stieltjes quasi-polynomials 

We consider the Fuchsian differential equation 



3=0 



k ,1 

a 3 ' 1 I 
1 V 



t — dj 
3=0 3 



+ 



1 

4e 



Pj a j A j 
t — dj 



k-l 



v = 



(2.1) 



for the function v(t) where the prime denotes differentiation with respect to t. This differential 
equation contains four sets of real parameters: 



a < at < ■ ■ ■ < ajfc, 
a ,ai, ...,a k e (-|, oo), 
Po,Pi> ■ ■ ■ iPk ^ {0, 1}, 
Aq, Ai, . . . Afc_i G M, 



and A,- is an abbreviation: 



A 3 := H(c 



(2.2) 
(2.3) 
(2.4) 
(2.5) 



(2.6) 



i=0 



Usually, the first three sets of parameters are given while the A's play the role of eigenvalue 
parameters. 

Equation (|2.1|) has regular singularities at infinity and at each ctj, j = 0,1,..., fc. The 
exponents at aj are Vj = and fij = —Jr 1 - — aj. If Uk+i, fJ-k+l denote the exponents at infinity 
then 



and 



Afe-l — Vk+lUk+l 



fc+1 

J2(uj+fij) = k. 

j=0 



(2.7) 



(2i 



The accessory parameters Ao, Ai, . . . , Xk-2 are unrelated to the exponents. 

The following result defines Stieltjes quasi-polynomials E n ^ p . Let parameter sets (12.21) and 
(|2,3|) be given. For every multi-index n = (ni, rt2, • • • , n^) G N fc of nonnegative integers and p = 
(PO)Pii • • • iPk) S {0, l} k there exist uniquely determined values of the parameters Ao, . . . , Afc_i 
such that (12.11) admits a solution of the form 



E n , P (t) 




\ p il 2 E n>p (t), t G M, 



(2.9) 



where E np is a polynomial with exactly rij zeros in the open interval (a,-_i,a,-) for each j = 
1, . . . , k. The polynomial E n>p (t) is uniquely determined up to a constant factor and has the 
degree \n\ = n± + ■ ■ ■ + n&. We normalize E np so that its leading coefficient is unity. Then we 
may write E nj) in the form 



E n , P (t) = Y[\t-a j \*/ 2 Y[{t-e t 

j=0 1=1 



(2.10) 
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where 

01 <9 2 < ■■■ < 0\ n \- 

Then 0\, . . . 9 ni lie in (a , a{), 9 ni +i, ■ ■ ■ ,0 ni +n 2 lie in (01, o 2 ) and so on. 

If p = then -E^o is a polynomial introduced by Stieltjes whose work was influenced 
by Heine Part III]. A proof of existence and uniqueness of the polynomials E n $ can be 
found in Szego ^3 Section 6.8]. For general p a computation shows that E ntP is the Stieltjes 
polynomial E n Q with ctj replaced by ay + pj. Therefore, the proof of existence and uniqueness 
in the general case can be reduced to the special case p = 0. 

The value of Xk-i associated with E np can be computed. One of the exponents at infinity 
must be 

k 

u k+1 = \n\ + \p\:=^2pj. 

j=o 

Using 11277)1 . (|275)l . we obtain 

k 

Afe_i = — \m (|m + |ce| + , m := 2|n| + |a| := ^Jay. (2-11) 

i=o 

No formulas are known for the corresponding values of the accessory parameters. 

If Qj = for all j then Stieltjes quasi-polynomials reduce to Lame quasi-polynomials in 
arbitrary dimension. If we work in R 3 there are eight possible choices of the parameters (|2.4|) 
giving us the familiar eight types of classical Lame quasi-polynomials; see Arscott pQ. 

3 Generalized ellipsoidal harmonics 

We say that a function u : R fc+1 — > R has parity p = (po, . . . ,p k ) £ {0, l} fe+1 if 
u{x) — u(ojx) = 2pju(x) for j = 0, 1, . . . , k. 
Equation l|1.4[) can be written in the form 

k „ rs k 



Au{x) + ^-^-u{x) - £ ^(«(a;) - ufosc)) = 0. (3.1) 

j=0 3 3 j=0 3 

If u has parity p then (j3.1[) becomes the partial differential equation 

Au( X ) + 2 ^ U{X) ~ £ = °- (3 ' 2) 

j=0 3 3 j=0 3 

In order to introduce ellipsoidal coordinates, fix the parameters (|2.2|) . For every (xo, . ■ ■ , xj-) 
in the positive cone of R fc+1 

x > 0, . . . ,Xk > 0, (3.3) 

its ellipsoidal coordinates to, t\, . . . , t k lie in the intervals 

a k < t < oo, cii-i <ti < ai, i = l,...,k, (3.4) 
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and satisfy 

k x 2 

— = 1 for i = 0,l,...,k. (3.5) 

Conversely, for given ti in the intervals (j3.4j) . we have 

k 

x) = *f . (3.6) 

ft (fk - aj) 

i=0 

These coordinates provide a bijective mapping between the positive cone (|3~3jl and the cube lEOIl . 

We now transform the partial differential equation (|3.2|) for functions u(x) defined on the 
cone (|3.3|l to ellipsoidal coordinates, and then we apply the method of separation of vari- 
ables. We obtain k + 1 times the Fuchsian equation (|2.1j) coupled by the separation constants 
Ao, Ai, . . . , Afc_i. We do not carry out the details of these known calculations. A good reference 
is Schmidt and Wolf Therefore, if Vj(tj), j = 0,1..., k, are solutions of (|2.1|) with tj 

ranging in the intervals (|3.4|) then the function 

u(x , ..., Xk ) = v (t ) ■ --Vkitk) (3.7) 

satisfies ()3.2|) . Of course, the values of the parameter sets (|2.2jl - (|2.5jl must be the same in each 
equation (|2.1|) . 

As a special case choose Vj as the Stieltjes quasi-polynomial E n p for each j. Then we know 
that 

xi,...,x k ) := E n>p (t )E n>p (ti) ■ ■ ■ En )P (t k ) (3.8) 
solves (|3.2|) . This function F njP is our generalized ellipsoidal harmonic. 

Theorem 1. The generalized ellipsoidal harmonic F np is a polynomial in xq,x\, . . . ,x k which 
satisfies Dunkl's equation (|1.4j) . It is of total degree 2\n\ + |p| and has parity p. 

Proof. If to, . . . , t k denote ellipsoidal coordinates of xq, . . . ,x k , then 

k / k \ ( k 2 

\{{h-o)= [Xl^-e)) i-E^rl ( 3 - 9 ) 

j=0 \i=0 / \ j=0 

for every 9 different from each aj. In fact, both sides of ()3.9|) are polynomials in 6 of degree k + 1 
with leading coefficient (— Moreover, both sides of the equation vanish at 9 = to, t\, . . . ,t k 
by definition ()3.5|) . So equation ()3.9|) follows. 

By (PHU)) . (j3~o1) . (|^H|) and (j33|) . the function F np can be written as 




F n , P (x) = c n<p x p I I > — ^ 1 , (3.10) 



where 

x p := xl° ■ ■ ■ x p k k for x = (x , . . . ,x k ), p = (p , . . . ,p k ), 
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and c nt p is the constant 

^ P :=(-i) w (ni^-r 72 ) nn^-^) (3.H) 

\j=o J \e=ii=o J 

with Aj according to (|2.6|) . This shows that F np is a polynomial of total degree 2\n\ + We 
know that F np solves ()3.2)) on the cone ()3.3)) and since it has parity p it solves (|1.4|) on ■ 



4 Generalized sphero-conal harmonics 

In order to introduce sphero-conal coordinates, fix the parameters (|2.2|) . Let (xq, xi, . . . , Xk) be 
in the positive cone of R fc+1 . Its sphero-conal coordinates r, s\, . . . , s k are determined in 
the intervals 

r > 0, aj_i < Sj < z = 1, . . . , k (4-1) 
by the equations 

r 2 = X>f (4.2) 

j=0 



and 



fe 2 



— = for i = 1, . . . , k. (4.3) 

j=o Si ~ aj 

This defines a bijective map from the positive cone in to the set of points (r, s\, . . . , s k ) 

satisfying ()4.1() . The inverse map is given by 

k 

x)=r^ . (4.4) 

n («i - a i) 

We now transforni th.6 pctrticil differentictl equation (j3.2|) for functions ii(xo ; . . . , ) defined 
on the cone ()3.3)) to sphero-conal coordinates and then we apply the method of separation of 
variables |TJS]. For the variable r we obtain the Euler equation 

// . k + 2\a\ , 4A fc _i 
u H v H ^—^0 = (4.5) 

while for the variables s\,S2, ■ ■ ■ ,s k we obtain the Fuchsian equation (|2.1|) . More precisely, if 
Ao, • • • , Afc_i are any given numbers (separation constants), if vo(r), r > 0, solves 1)4. 5 |l and Vi(si), 
Oj_i < Sj < <2j, solve 1)2.1)) for each i = 1, . . . , k, then 

n(x ,a;i, ...,x k )= v (r)vi(si)v 2 (s 2 ) ■ • -ffc(sfc) 
solves (EU). 

Let -En iP be a Stieltjes quasi-polynomial. It follows from 1)2.111) that uo(?") = r m is a solution 
of 1)4. 5|) . where m := 2|n| + |p|. Therefore, 

G^o.ii,...,^) := r m E n)P (s 1 )E ri) p(s2) ■ ■ ■ E n>p (s k ) (4.6) 
is a solution of ()3.2)) . This function G TCiP is our generalized sphero-conal harmonic. 
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Theorem 2. The generalized sphero-conal harmonic G n ^ p is a polynomial in xq,x\, . . . ,x k , it 
is homogeneous of degree 2\n\ + \p\, it has parity p and it solves Dunkl's equation (|1.4|) . 

Proof. Let (xq, . . . , x k ) be a point with xj > for all j, and let (r, s\, . . . , s^) denote its 
corresponding sphero-conal coordinates. We claim that 



r\ Sl -8)...(s k -0)=[ TTK - 9) V —L- (4.7) 




for all 8 which are different from each aj . Both sides of (|4.7|) are polynomials in of degree k 
with leading coefficient (— l) k r 2 . Moreover, both sides vanish at 6 = si,...,sj. because of 
definition (|4.3j) . Equation (|4.7j) is established. 

We write the Stieltjes quasi-polynomial E np in the form (j2. 1U|) . Using (|4.4f) . ()4.6|) and l|4.7|) . 
we obtain 

I "I fe ^.2 

G„, p (a:) = TT V , (4.8) 

where c ntP is given by (|3.11j) . This shows that G ntP (x) is a polynomial in xo,xi, . . . , x^, it is 
homogeneous of degree 2\n\ + and it has parity p. We know that G ntP solves ()3.2(l on the 
cone ()3.3(l and since it has parity p it solves (|1.4j) on R fc+1 . ■ 

A generalized spherical harmonic is a homogeneous polynomial u in the variables xq, x±, . . . , x\~ 
which solves Dunkl's equation (|1.4|) . For a given set of parameters (|2.3|) we let 7i m denote the 
finite dimensional linear space of all generalized spherical harmonics of degree m. If aj = for 
each j then we obtain the classical spherical harmonics. 

On the fe-dimensional unit sphere S k we introduce the inner product 



(f,g) w -= w(x)f(x)g(x)dS(x), (4.9) 
Js k 

and norm 

ll/L:=(/,/)l /2 , (4-10) 

where the weight function w is defined by 

w(x , x u ...,x k ):= |x | 2Q0 N 2ai • • • \x k \ 2a K (4.11) 

The surface measure on the sphere is normalized so that J S k dS(x) equals the surface area of the 
sphere S k . The condition ay > — ^ ensures that {f,g) w is well-defined if / and g are continuous 
on S k . 

Theorem 3. Let m G N. The system of all generalized sphero-conal harmonics G n ^ p of degree m 
forms an orthogonal basis for 7i m with respect to the inner product ()4.9|) . 

Proof. We consider the system of all sphero-conal harmonics G ntP , where n,p satisfy m = 
2\n\ + \p\. By Theorem|2j G ntP belongs to 7i m . The dimension of the linear space of generalized 
spherical harmonics of degree m which have parity p is 



- 2 |p| + k — 1 
k - 1 



(4.12) 



if m — \p\ is a nonnegative even integer and zero otherwise. This can be proved as in Hochstadt 
[71 p. 170] or it follows from Dunkl [31 Proposition 2.6] where a basis of TL m in terms of Jacobi 
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polynomials is constructed. The dimension Q4.12JI agrees with the number of multi-indices 
n = (m, ri2, ■ ■ ■ , Tift) for which m = 2\n\ + \p\. We conclude that the number of pairs n,p 
with m = 2|n| + \p\ agrees with the dimension of 7i m . Therefore, in order to complete the proof 
of the theorem, we have to show that G nyP is orthogonal to G n tw provided (n,p) ^ (n',p'). If 
p ^ p' this is clear because the weight function (|4.11|) is an even function. If p = and n ^ n' 
then orthogonality was shown in |19| Theorem 3.3]. The proof of orthogonality in the remaining 
cases is analogous and is omitted. ■ 

Extending the method of proof of |19| Theorem 3.3] we also establish the following theorem. 

Theorem 4. The system of all generalized sphero-conal harmonics G n ^ p , n 6 N fe , p € {0, 
when properly normalized, forms an orthonormal basis of L^ u (S k ). 

No explicit formula is known for the norm of G np in Lf IJ (S k ). However, the norm of a poly- 
nomial can be computed using the formula 

2 n m) 

Nr-vr- 1 • • • w*- 1 ds(x) = — - (4.i3) 

s k r (Po + Pi H h Pk) 

which holds whenever f3j > 0, j = 0, 1, . . . , k. 

In general, ellipsoidal harmonics are not homogeneous polynomials so they are not spherical 
harmonics. However, they are related to spherical harmonics in the following way. 

Theorem 5. Let m G N. The system of all generalized ellipsoidal harmonics F ntP of total degree 
2\n\ + \p\ at most m is a basis for the direct sum 

H ®Hi®---®H m . (4.14) 

Proof. The Dunkl operator A a maps a homogeneous polynomial of degree q to a homogeneous 
polynomial of degree q — 2. Therefore, if we write a generalized ellipsoidal harmonic as a sum of 
homogeneous polynomials, then these homogeneous polynomials also satisfy (|1.4j) . Hence every 
generalized ellipsoidal harmonic of total degree at most m lies in the direct sum (|4.14|) . By 
comparing l|3.1Ujl with (|4.8|) . we see that 

F n , P (x) = G n)P (x) + terms of lower degree. (4-15) 

By Theorem |3 the system of all G niP with 2|n| + |p| < m is a basis for the direct sum (|4.14|) . 
The statement of the theorem follows. ■ 

Of course, the spaces 7i m of generalized spherical harmonics and the direct sum (|4.14j) depend 
on the parameters ao,a\, . . . , a^. However, it is easy to show that the set of functions on S k 
which are restrictions of function in (|4.14|) is independent of these parameters. In fact, this set 
consists of all functions that are restrictions of polynomials of total degree at most m to S k . 

5 Hobson's formulas 

In this section we generalize some formulas given by Hobson |HJ p. 124]. These formulas will be 
applied in the next section to obtain a generalization of Niven's formula. 

Lemma 1. Let T> be the operator given by 
Vf{x) := f (x) + a , 
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where a is a constant. Then, for A(z) := z , f,m6N, we have 

m ^ 

v m x 2i = 2 m -%A<- m ~fi (x 2 ) -V 2j x m . (5.1) 

3=0 J - 

Proof. If m > 2£, then both sides of 1)5.1)) are zero. So we assume that m < 2£. We first 
consider the case that m = 2n is even. The left-hand side of (|5.1j) is equal to 

2 m (^(-ir(^-<^" m (5-2) 

The right-hand side of (|5,lj) is equal to 

Jl 1 2 2,- n\ 

m + j)\j\ (n-j)\ 



it 



y 2 m ~^—^——-2^ T ^—-{-\y a-a-n). x 2£ ~ m . (5.3) 

j=0 



After some simplifications, equality of ()5.2)1 and (|5.3f) follows from the Chu-Vandermonde sum 

A (a)j(b)n-j _ (a + fe)n 

applied to a = ^ — a — n, b = n — £. This completes the proof of 1)5.1)) if m is even. The similar 
proof for odd m is omitted. ■ 

Clearly, in 1)5.1(1 it would be enough to let j run from to [^J • Similar remarks apply to 
other formulas in this section. 

In the following lemma, T>j is according to 1)1.5)1 and dj is the usual partial derivative with 
respect to Xj. 

Lemma 2. Let mo, mi, . . . , m& 6 N, and Ze£ A : (0, oo) fc+1 R k ra := mo + mi + • • • + m^ 
times differentiable. Then, for xq, . . . ,Xk ^ 0, 

mo jnk 

io=o j' fe =o 

x (^°-^ . . . fl^-'MKxg, . . . , 4) p v" p r ^° • • • < fe - ( 5 - 4 ) 

Jo! • • • Jfc! 

Warning: On the left-hand side of this formula the operators T>j are applied to the function 
f(xo,...,Xk) '■= A(xq, . . . ,x 2 ), whereas on the right-hand side the partial derivatives dj are 
applied directly to A. 

Proof. Let B be the Taylor polynomial of A of order m at a given point {zq, . . . , z^) with Zj > 0. 

1 /2 M 

Let Xj := Zj . Then 1)5.4)1 is true with ^4 — B in place of A at the point Xo, ■ ■ ■ ,x k (both sides 
of the equation are zero.) Therefore, it is sufficient to prove 1)5. 4|) for polynomials A, and so for 
monomials 



A(z Q ,...,z k ) = z 



-k ■ 



In this case, we obtain 1)5.4)1 by applying Lemma ^ to each function z^ j and multiplying. ■ 

If /(xq, . . . , Xf~) is a polynomial, we will use the operator /(X>o> ■ ■ ■ >2?fc)- It is well-defined 
because the operators T>j commute. We use 



r := (x Q + ...+x k ) 



2n1/2 
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Theorem 6. Let f m (xo, . . . , x k ) be a homogeneous polynomial of degree m, and let B : (0, oo) — > 
R be m times differentiable. Then, for all nonzero (xo, x\, . . . ,x k ), 

f m (V , V k )[B(r 2 )) = 2 m -^B^\r 2 )-AU m (x , x k ). (5.5) 

3=0 J ' 

Proof. It is sufficient to prove (|5.5|) for monomials 

f m (x , ...,x k ) = x™° ■ ■■x k nk ) m = m -\ h m k . 

In this case (|5.5|) follows from Lemma|2]with A(zq, . . . , z k ) = B(zq + • • • + z k ) by using 

grno-jo . . . fl^-^fo, . . . , Z k ) = B^\z + ••• + «*) 
with j = jo + • • • + j k and the multinomial formula 



■ , , ■ _ • \J0 ' ' ' Jk/ 



joH hjfc=J 

When we apply Theorem H3 to B(z) := z 1 with 7 € K, we obtain the following corollary. 
Corollary 1. Let f m (xo, ■ ■ ■ ,x k ) be a homogeneous polynomial of degree m. Then, for 7 € M, 

m 

f m (V Q , . . . ,V k )[r 2 ^} = ^ 2 -2 J( _ ir - J( _ 7)m _ jr 2( 7 - m+J )_ AJJm(l()i (5 6) 

3=0 ' 7 ' 

Corollary 2. Let f m be a generalized spherical harmonic of degree m. Then, for 7 £ M, 

r 2(m " 7) f m (Po, ■ ■ .,V k )[r 2 ^] = 2 m (-l) m (- 7 ) m / m (xo, . . . , x k ). (5.7) 
If we set 

7 = 1 T ^-I«|, (5.8) 

then a simple calculation shows that A Q r 27 = 0. So r 27 plays the role of a fundamental solution 
of A a u = generalizing the solution 1/r of the Laplace equation Au = in M 3 with a = (0, 0, 0). 
Note that the number 7 defined by ()5.8|) is always less than 1. It can be zero (for example for 
the Laplacian in the plane). In this case, lnr plays the role of a fundamental solution. The 
fundamental solution r 27 and an associated formula producing harmonic polynomials appeared 
in Xu E2. 



Corollary 3. Lf / m (xo, • • • ,x k ) is a homogeneous polynomial of degree m and 7 is defined by 
t)5.8|) . then the right-hand side of equation (|5.6j) is a generalized spherical harmonic of degree m. 

Proof. Since A a r 21 = 0, this follows by applying A a to both sides of (|5.fi[) . ■ 

6 Niven's formula 

In this section we prove a generalization of Niven's formula expressing ellipsoidal harmonics in 
terms of sphero-conal harmonics. We follow the method of Hobson p. 483]. 
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Let E = E n ^ p be a Stieltjes quasi-polynomial, and let F = F np , G = G njP be the correspond- 
ing ellipsoidal and sphero-conal harmonics written in the forms (|3.1U|) and (|4.8|) . respectively. It 
will be convenient to introduce the auxiliary polynomial 



x 2 Ji, x 2 



H(x) :=c n , pX PH E^TT-E 



t — cij 

where t is a fixed number greater than a k . We define positive constants dj by d 2 = t — aj for 
j = 0, . . . , k. Let 7 be the constant defined by (|5.8j) . We assume that 7 / 0. The identity 

k x 2 - k x 2 k x 2 

E I^a~ ~ E T^a~ = (t ~ ^ E (Q-aAt-ajY 

j=0 J j=0 J j=Q v J ' x Jl 

implies 

H(d x , d k x k ) = E(t)G(x , x k ). (6.1) 
By Corollary [21 we have 

2 m (-l) m (- 1 ) m G(x , ...,x k ) = r 2 ^G(V , . . . ,V k )[r 2 ^ (6.2) 
where m := 2\n\ + \p\. Since 

k x 2 _ k (t - aj )x 2 



r 2 



^-*)E^r = E 0-a 3 

j=Q J j=Q J 



we conclude that 

G(d x , d k x k ) = E(t)G(x , ...,x k )+ r 2 P(x , x k ), 
where P is a polynomial. It follows that 

G(d V , d k V k ) = E(t)G(V , ...,V k )+ P{V Q , V k )A a . 
Using that A a r 27 = 0, we obtain 

G(d V , d k V k ) [r 2 ^] = E(t)G(V , V k ) [r 2 ^] . (6.3) 
We now combine equations (|6.1|) . (|6.2|) . ()6.3|) and obtain 

2 m (-l) m (- 7 ) m H(d x , . . .,d k x k ) = r 2 ^G(d V , . . . ,d k V k )[r^] . (6.4) 
Now (|6.4|) and Corollary ^ yield 

™ r 2i 

H(d x , ■ ■ -,d k x k ) = } —-— — — A l a [G(d Q xo, ■ ■ -,d k x k )}. 

f^J 2 z Hl{j — m+ l)i 

We replace the variables Xj by Uj/dj and rename yj as Xj again. This gives 

m r>2i 

H W = E 22 ^ !(7 _ m + 1 v (^o 2 + • • • + 4n) l G(x), (6.5) 

i=0 



where 



k 2 k 2 



* 2 = £3 = £ 



Since G satisfies A Q G = 0, we see that the right-hand side of (|6.5|) does not change if we replace 
d 2 V 2 by —ajV 2 . If positive numbers xq,... ,x k are given, we can choose t as the ellipsoidal 
coordinate t = to. Then R = 1 and H(x) = F{x). We have proved the following theorem. 
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Theorem 7. Let F n ^ p , G n ^ p be the generalized ellipsoidal and sphero-conal harmonics of degree 
m = 2\n\ + \p\ defined by (|3.1U|) . (|4.8|) . respectively. Assume that 7 defined by (|5.8jl is nonzero. 
Then 

m (—lY 

F ^ x ) = E 2M7 + l - m) > ° P o + • • • + a k VlYG n , p (x). (6.6) 



i=0 



In the classical case k = 2 and oeo = a\ = 02 = this is Niven's formula; see [HJ p. 489]. 



7 A Dirichlet problem for ellipsoids 

In this section we apply generalized ellipsoidal harmonics to solve the Dirichlet boundary value 
problem for the Dunkl equation on ellipsoids. 
We consider the solid ellipsoid 



£ :-- 



k x 2 

x £ R k+1 : ^ -i < 1 
3=0 b i 



with semi-axes bo > b± > ■ ■ ■ > b k > 0. Let d£ be the boundary of £. Given a function 
/ : d£ — > M we want to find a solution u of Dunkl's equation ljl.4J) on £ that assumes the given 
boundary values / on d£ in the sense explained below. 

It will be convenient to parameterize d£ by the unit sphere S k employing the map 

T :S k -> d£ (7.1) 

defined by 

T(yo,yi, ■ ■ ■ ,Vk) ■= {boyo,bwi,...,b k y k ). 

We suppose that the given boundary value function / : d£ — ► M has the property that the 
function / o T is in L 2 v (S k ), where the weight function w is defined in ()4.11() . A solution of the 
Dirichlet boundary value problem for the Dunkl equation with given boundary value function / 
is a function u £ C 2 (£) which satisfies (|1.4[l in £ and assumes the boundary value / in the 
following sense. For sufficiently small 5 > 0, form the confocal ellipsoids 

*^ k+i --i: A=4> (7.2) 
3=0 3 ) 

and let T$ be defined as T but with respect to the ellipsoid (|7.2|) in place of d£ . Then we require 
that 

uoTs^foT in L 2 w (S k ) as 5^0. (7.3) 

We now show how to construct a solution of this Dirichlet problem. We choose any real 
number uj (we can take uj = if we wish), and define numbers ao < a\ < ■ ■ ■ < by 

aj := uj — b 2 . 

Corresponding to these numbers aj we introduce sphero-conal coordinates (r, si, . . . , s k ) for 
cartesian coordinates y = (yo, y\, . . . , y k ) and ellipsoidal coordinates (to, <i, . . . , tfc) for cartesian 
coordinates x = (xq, x±, . . . , x^j. Note that if x = Ty and r = 1 then Sj = tj for each 
j = l,2,...,k. 
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Since the function / o T lies in L 2 L ,(S k ), we can expand / o T in the orthonormal basis of 
Theorem |1J 

f oT = ^2 fn,P e n, P Gn,p, (7.4) 

71, p 

where the factors e n ^ p are determined by 
^n,p||Cn,p||iu 1, 

and 

fn,p • — (/ ° T) 6n,pGn t p) w (7.5) 

Then 

^|/„, p | 2 = ||/oT||2 <oo. (7.6) 

n,p 

The expansion (|7.4j) converges in L 2 u (S k ). We are going to prove that 

u(x) := £ (a;) (7 . 7 ) 

is the desired solution of our Dirichlet problem. 

Theorem 8. The function u defined by (|7,7|) is infinitely many times differentiate and solves 
Dunkl's equation (|1.4j) on the open ellipsoid £ , and it assumes the given boundary value f in the 
sense of (|7.3|) . 

Proof. We first show that the infinite series in l|7.7|) converges. We know from (141 Lemma 2.2] 
that there is a sequence K m of polynomial growth such that 

\e n , P G n ,p{y)\ < K m for m = 2\n\ + \p\, y £ S k . (7.8) 
For given t £ {ak,uj) we consider the solid ellipsoid 



£ 



t ■- 



k X 2 

x £ R k+1 : V J — < 1 

^ t-a-i 



which is a subset of £. By comparing (|3.8|) and l|4.6|) . we get from (|7.8|) 

\z n ,pF n , p {x)\ < E njP (t)K m for x £ d£ t . (7.9) 

The Stieltjes quasi-polynomial E n p has degree m/2 and all of its zeros lie in the interval [ao, af-]- 
Hence we have the inequality 

— — ^ ) for a k < t < uj. (7.10) 

u — ao / 

Now we obtain from (f7T9|) . (fTTTUl) 

(t a \ m / 2 
) E ntP (uj)K m for x £ £ t . (7.11) 
w — ao/ 
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Since the set of numbers f niP is bounded by (|7,6|) . K m grows only polynomially with m and 

( ili-ao ) m ^ S oes to exponentially as m — ► oo, we see that the series in (|7.7|) converges uniformly 
in £ t and thus in every compact subset of £. 

The next step is to show that u solves equation (|1.4|) . This follows if we can justify inter- 
changing the operator A a with the sum in (|7.7j) . Consider the series 

that we obtain from (|7.7|) by differentiating each term with respect to Xj. In order to show 
uniform convergence of this series on St we need a bound for the partial derivatives of F np . 
We obtain such a bound from the following result due to Kellog If P(yo, Hi, ■ ■ ■ , Uk) is a 
polynomial of total degree N then 

max{||gradP(y)|| : \\y\\ < 1} < iV 2 max{|P(y)| : \\y\\ < 1}, (7.13) 

where || • || denotes euclidian norm in If we use a mapping like (|7.1j) to transform the 

ellipsoid to the unit ball we find that 



d 

Fn )P \X) 



d Xj 



9 

m 



^ n ^ max{|F T , jP (z)| : z G S t ] for x G £ t . (7.14) 



Using this estimate we show as before that ()7.12(l converges uniformly on £f In a similar way 
we argue for the second term in the generalized partial derivative (jl.5j) . Since we can repeat 
the procedure we see that u is infinitely many times differentiable on £ and it solves Dunkl's 
equation. 

It remains to show that u satisfies the boundary condition (|7.3|) . For y £ S k , we find 
/ o T{y) - f o T s (y) = £ f n , p e n , p (l - E ^ ~ S A G ( y ). 

Hence we obtain 

ii/ o t - f o Ts \\i = y: n, P (i - 2 

so ()7.3|) follows easily. ■ 

8 Examples 

Formulas in this paper have been checked with the software Maple for some Stieltjes polynomials 
represented in explicit form. For example, take k = 2, 

229 71 25 

a = 0, ai = 3, a 2 = 5, a = — , «i = 77, «2 = 

54 54 

and 

ni = 2, n 2 = 1, Po = Pi=P2 = 0. 
Then the corresponding Stieltjes polynomial E n ^ p is given by 
E n>p (t) = (t-l)(t-2)(t-4). 
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Indeed, E np satisfies equation (|2.1j) with 

1120 119 

Ao - — , Ai = - — , 

and it has two zeros between ao and oi, and one zero between a\ and 02- 

The simplest way to compute such examples is to use the fact that the zeros 0£ of E n o are 
characterized by the system of equations 



\ -» 2 \ a j + 2 



9=1 * 9 7=0 * 3 

see (H (6.81.5)]. 

The corresponding ellipsoidal and sphero-conal harmonics are 

771 _1 Q9 / 2 _ 1 2 1 2 -i\ (1 2 2 1 2 i\ (1 2 1 2 2 -A 

^n,p — i-V^ [Xq 2^1 4 ;E 2 X J i.^O X l 3 2 '2 X J \4 X ' X l X 2 L ) > 

/° TOO/' 2 l 2 12W12 2 12W12.2 2\ 

^n.p — — — 2 X 1 ~~ ~i x 2) \2 X ~ x l ~ 3 X 2) U^O + x l ~ X 2J ■ 

One can check that these polynomials do satisfy equation (|1.4|) . Also, applying formula (|6.6f) to 
G„ iP we obtain F n p as claimed. 

We now take the same aj but replace the parameters aj by 

229 175 71 17 25 

00 = "54 " 1 = "54 ' ai = 54 - 1 = 54' ° 3 = y 

Moreover, let 

ni = 2, ra 2 = 1, po = 1, Pi = 1, P2 = 0. 
Then the Stieltjes quasi-polynomial is 



£n, P (t) = VI*IVI* - - l)(t - 2)(t - 4). 
It satisfies equation (|2,lj) with 

2855 440 

The corresponding ellipsoidal and sphero-conal harmonics are as before but with —192 replaced 
by — 192v / 15v / 6 and the extra factor x±X2 added. Again it can be checked that these polynomials 
satisfy equation Q1.4JI . and formula 1)6.6(1 holds. 
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